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Secular increase of the Astronomical Unit: 
a possible explanation in terms of 
the total angular momentum conservation law 
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Abstract 

We give an idea and the order-of-magnitude estimations to explain the recently reported secular increase 
of the Astronomical Unit (AU) by Krasinsky and Brumberg (2004). The idea proposed is analogous to 
the tidal acceleration in the Earth-Moon system, which is based on the conservation of the total angular 
momentum and we apply this scenario to the Sun-planets system. Assuming the existence of some tidal 
interactions that transfer the rotational angular momentum of the Sun and using reported value of the 
positive secular trend in the astronomical unit, ^AU = 15±4(m/cy), the suggested change in the period 
of rotation of the Sun is about 21ms/cy in the case that the orbits of the eight planets have the same 
"expansion rate." This value is sufficiently small, and at present it seems there are no observational data 
which exclude this possibility. Effects of the change in the Sun's moment of inertia is also investigated. It 
is pointed out that the change in the moment of inertia due to the radiative mass loss by the Sun may 
be responsible for the secular increase of AU, if the orbital "expansion" is happening only in the inner 
planets system. Although the existence of some tidal interactions is assumed between the Sun and planets, 
concrete mechanisms of the angular momentum transfer are not discussed in this paper, which remain to 
be done as future investigations. 
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1. Introduction 

The Astronomical Unit (hereafter we abbreviate AU) is 
one of the most essential scale in astronomy which char- 
acterizes the scale of the solar system and the standard of 
cosmological distance ladder. AU is also the fundamental 
astronomical constant that associates two length unit; one 
(m) in International System (SI) of Units and one (AU) 
in Astronomical System of Units. 

In the field of fundamental astronomy e.g., the plane- 
tary ephemerides, it is one of the most important subjects 
to evaluate AU from the observational data. The appear- 
ance of planetary radar and spacecraft ranging techniques 
has led to great improvement of determination of AU. 
Modern observations of the major planets, including the 
planetary exploration spacecrafts such as Martian landers 
and orbiters, make it possible to control the value of AU 
within the accuracy of one meter or even better. Actually 
present best-fit value of AU is obtained as (Pitjeva 2005), 

1 (AU) = AU (m)== 1.495978706960 xl0"±0.1 (m).(l) 

However, recently Krasinsky and Brumberg reported 
the positive secular trend in AU as ^AU = 15±4(m/ey) 
^ from the analysis of radar ranging of inner planets and 
Martian landers and orbiters (Krasinsky and Brumberg 



In this paper, cy means the century according 
Brumberg (Krasinsky and Brumberg 2004). 



to Krasinsky and 



2004; Standish 2005). As mentioned above, the value of 
AU is presently obtained within the error of one meter or 
even better, so that the reported dAU/dt — 15±4(m/cy) 
seems to be in excess of current determination error of 
AU. This secular increase of AU was discovered by using 
following formula (Krasinsky 2007), 
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AU + -— i-to 
dt 
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in which ^theo is the theoretical value of round-trip time 
of radar signal in the SI second, dtheo is interplanetary 
distance obtained from ephemerides in the unit of (AU), 
c is the speed of light in vacuum, AU and dASJ / dt are, 
respectively, the astronomical unit as Eq. (1) and its time 
variation, and t^ is the initial epoch of ephemerides. tthco 
is compared with the observed lapse time of signal <obs 
and then AU, dAU/dt are fit by the least square method. 

It should be emphasized that, at this time, the exis- 
tence of secular increase of AU is not confirmed robustly 
by other ephemerides groups, and there is still room for 
the independent analysis of observational data. However, 
several attempts have already been made to explain this 
possible secular increase of AU, including e.g., the ef- 
fects of the cosmic expansion (Krasinsky and Brumberg 
2004; Mashhoon et al. 2007; Arakida 2009), mass loss 
of the Sun (Krasinsky and Brumberg 2004, Noerdlinger 
2008), the time variation of gravitational constant G 
(Krasinsky and Brumberg 2004), the influence of dark 
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matter (Arakida 2008) and so on. But unfortunately so 
far none of them seems to be successful. 

In this paper, we will take another viewpoint and give 
an explanation for the secular increase of AU, based on 
the standard conservation law of the total angular mo- 
mentum. The mechanism is closely analogous to the case 
of the tidal acceleration in the Earth-Moon system which 
induces the transfer of the rotational angular momentum 
of the Earth into the orbital angular momentum of the 
Moon due to tidal friction. We will apply similar scenario 
to the Sun-planets system. 

2. Tidal acceleration in the Earth-Moon system: 
a brief summary 

In this section, we briefly summarize the tidal acceler- 
ation in the Earth-Moon system. See e.g.. Chapter 4 of 
Murray and Dermott (1999) for details. 

The conservation law of the total angular momentum 
in the Earth-Moon system is 

^(£E + iM) = 0, (3) 

where £e is the rotational angular momentum of the 
Earth, Lm is the orbital angular momentum of the Moon, 
and for simplicity we have neglected the rotational angu- 
lar momentum of the Moon, which is about 10~^ times 
smaller than Lm- The rotational angular momentum is 
written in terms of the moment of inertia and the pe- 
riod of rotation Te as 

^e = /e^. (4) 
The moment of inertia can be written as 



(5) 



where M is the mass and R is the radius of the Earth, 
and 7 is the moment of inertia factor. In the case of the 
perfect sphere with uniform density inside, 7 is 2/5. For 
a centrally-condensed body, 7 is in general less than 2/5. 
Often quoted value for the Earth is 7 = 0.3308. The orbital 
angular momentum of the Moon is written as 



Lm = mVGA/r(l-e2), 



(6) 



where m is the mass of the Moon, r and e is the semi- 
major axis and the eccentricity of the Moon's orbit, re- 
spectively, and G is the gravitational constant. For the 
sake of simplicity, hereafter we assume e is constant. 

Then, if we assume /e, m, and M are constant, the 
conservation of the total angular momentum Eq. (3) reads 



(7) 



Te _ 1 Lm r 
Te " 2 ^E r 

The motion of the Moon can be traced with an accu- 
racy of a few centimeters by the lunar laser ranging ex- 
periments. The measurements yield the numerical value 
f = -f3.82 ± 0.07 (m/cy) (Dickey et al. 1994) then 



1.0 X 10-*(cy-i). 



(8) 



Using Eq. (8), the numerical value of Eq. (7) is 



Te 



.Ix 7-1 X 10-9 (cy-i). 



(9) 



or equivalently, Te — 2.1 (ms/cy), where we have used 7 = 
0.3308. 

On the other hand, historical records over 2700 years 
can be used to compute the change in the length of the 
day, and the following average value is found: 



-1.70 ±0.05 (ms/cy). 



(10) 



The gradual slowing of the Earth's rotation is due to the 
tidal force between the orbiting Moon and the Earth, or 
the tidal friction. It is explained that the discrepancy 
between the estimated value Te and the observed one Tobs 
is largely due to the change in the moment of inertia of 
the Earth caused by the melting of ice at the poles. 

3. AppHcation to the Sun-planets system 

In this section, we apply the same argument in the pre- 
vious section to the Sun-planets system. All we need is 
the conservation of the total angular momentum in the so- 
lar system. We denote the mass and the orbital elements 
of each planet by subscript i. The length AU, denoted 
by a, is used to normalize r^, then for the Earth's radius 
7-3 = Te = a, and for the moment it is assumed that the 
orbits of the all planets have the same "expansion rate," 



(11) 



rj _ a 
Ti a 
for each i. 

In the case that the Sun's moment of inertia does not 
change, from the analogy to Eq. (7), we obtain the change 
in the period of rotation of the Sun Tq as 

Tq 1 L a 

where L is the sum of the orbital angular momentums of 
all planets 



(12) 



,y'GMor,(l-ef), 



is the rotational angular momentum of the Sun 
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(14) 



and Af0, Rq, and 70 are the mass, the radius, and the 
moment of inertia factor of the Sun, respectively. 

Using the value d ~ 15(m/cy) reported by Krasinsky 
and Brumberg, d/a^ 1.0 x lO-^" and the right-hand-side 
of Eq. (12) is evaluated as 



p^^5.7x70^xlO-i°(cy-i). 



(15) 



If we use the value 7© = 0.0059 and the rotational period 
of the Sun as Tq = 25.38 (days), the estimated value of the 
change in Tq is 
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ro~21(ms/cy). 



(16) 5. Case of the Sun-inner planets system 



The estimated value is sufficiently small and seems to be 
well within the observational limit. 

4. EflFects of the change in the moment of inertia 

In the previous section, we have considered only the 
case that the moment of inertia Iq is constant, namely. 



dt 



d 
dt 



2tt 



IQ. 



(17) 



If we generalize our result to the case that /0 and Mq 
are not constant, Eq. (12) changes to 



70 _Mq _2i?Q T© 

70 Mq i?0 Tq 



1 L Mq a 
I Mq a 



(18) 



As a first approximation, we assume that the radiative 
mass loss occurs isotropically along radial direction and 
does not carry the angular momentum. 

The first term in the left-hand-side of Eq. (18) repre- 
sents the effect of change in the internal density distribu- 
tion of the Sun, and so far we do not have enough infor- 
mation on it in detail. 

The second term in the left-hand-side of Eq. (18) repre- 
sents the effect of mass loss, which can be evaluated in the 
following way. The Sun has luminosity at least 3.939 x 
10^^ W, or 4.382 x 10^ kg/s, including electromagnetic 
radiation and contribution from neutrinos (Noerdlinger 
2008). The particle mass loss rate by the solar wind is 
about 1.374 x 10^ kg/s, according to Noerdlinger (2008). 
The total solar mass loss rate is then 



= 9.1 X 10-12 (cy-l), 

Mq V ^ ^' 



(19) 



which is less than a thousandth of the required value 
to explain the secular increase of AU (see Eq. (15)). 
(As pointed out by Noerdlinger (2008), Krasinsky and 
Brumberg (2004) unaccountably ignored the radiative 
mass loss Lq = 3.86 x lO^^W which is the major contri- 
bution to Mq/Mq.) Therefore, we can conclude that the 
solar mass loss term in the left-hand-side of Eq. (18) does 
not make a significant contribution to the secular increase 
of AU, if the orbits of the eight planets have the same 
"expansion rate." 

Note that the term which is proportional to Mq/Mq 
also appears in the right-hand-side of Eq. (18). This 
term may be called as the Noerdlinger effect (Noerdlinger 
2008). Noerdlinger already investigated this effect of solar 
mass loss, and concluded that the effect can only account 
for less than a tenth of the reported value by Krasinsky 
and Brumberg. 

The third term in the left-hand-side of Eq. (18) is 
the contribution from the change in the solar radius. 
Although the very short-time and small-scale variability 
in the solar radius may be actually observed in the con- 
text of helioseismology, we have no detailed information 
on the secular change in Rq so far. 



In the previous sections, we have assumed that the or- 
bits of all the planets have the same "expansion rate", 
i.e., Eq. (11). However, the recent positional observa- 
tions of the planets with high accuracy are mostly done 
within the inner planets region. Actually Krasinsky and 
Brumberg obtained dAU/dt by using these inner planets 
data, while the orbits of outer planets are given by the 
Russian ephemeris EPM (Ephemerides of Planets and the 
Moon). Therefore, as a tentative approach, we consider 
the case that the "expansion rates" of the planetary or- 
bits are not homogeneous but inhomogeneous in the sense 
that only the orbits of the inner planets expand. 

In this section, we consider the case that the "expan- 
sion" of the planetary orbit occurs only for the inner plan- 
ets. 

In this case, the sum of the orbital angular momentum 
of all planets L is replaced by the sum of the inner planets 



Lin = ^ m, y'CAfe ( 1 - e| ) . 



Then Eq. (18) is now 

70 Mq ^Rq Tq Mq d" 

70 Mq Rq Tq 2 1q\Mq a 



(20) 



(21) 



Note that the sum of the angular momentum of inner 
planets amounts only 0.16% of the total L: 



Lin 



1.6 X 10" 



(22) 



Therefore, under the assumption that the change in the 
rotational angular momentum of the Sun affects only the 
orbital angular momentums of the inner planets, the re- 
quired values which were calculated in the previous sec- 
tions to explain the secular increase of AU can now be 
revised to be 1.6 x 10"^ times smaller. In particular, the 
right hand side of Eq. (21) is 



Jy^-- 1.5 X 10-11 (cy-i). 
2 Iq a 



(23) 



Interestingly, it is the same order of magnitude as (ac- 
tually it is about 1.6 times larger than) Mq/Mq (see 
Eq. (19)). 

Then we can conclude from Eq. (21) that the decrease 
of rotational angular momentum of the Sun due to the 
radiative mass loss has a significant contribution to the 
secular increase of the orbital radius of the inner planets. 

6. Conclusion 

In this paper, we considered the secular increase of 
astronomical unit recently reported by Krasinsky and 
Brumberg (2004), and suggested a possible explanation 
for this secular trend by means of the conservation law 
of total angular momentum. Assuming the existence of 
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some tidal interactions that transfer the angular momen- 
tum from the Sun to the planets system, we have obtained 
the following results. 

From the reported value ^AU= 15±4(m/cy), we have 
obtained the required value for the variation of rotational 
period of the Sun is about 21 (ms/cy), if we assume that 
eight planets in the solar system experience the same or- 
bital expansion rate. This value is sufficiently small, and 
at present it seems there are no observational data which 
exclude this possibility. 

Moreover, we have found that the effects of change in 
the moment of inertia of the Sun due to the radiative 
mass loss may be responsible for explaining the secular 
increase of AU. Especially, when we suppose that the or- 
bital expansion occurs only in the inner planets region, the 
decrease of rotational angular momentum of the Sun has 
enough contribution to the secular increase of the orbital 
radius. Then as an answer to the question "why is AU in- 
creasing ?" , we propose one possibility, namely "because 
the Sun is losing its angular momentum." 

As showed in previous section, the obtained value in 
Eq. (23) is practically about 1.6-times smaller than that 
due to Mq/Mq. Nonetheless the estimated dAU/dt in the 
different data sets and fitting parameters is distributed 
within the range 7.9±0.2 to 61.0±6.0 (m/cy) then seems 
not to be tightly constrained, see Table 2 of Krasinsky and 
Brumberg (2004). Hence we can say that the estimated 
value in Eq. (23) falls into the suitable range. 

In the process of planetary cphemeridcs construction, 
the effects due to both solar mass loss and rotational 
angular momentum transfer are currently not included. 
However, as mentioned by Noerdlinger (2008), the solar 
mass loss induces the variation in the orbital radius of 
planet i5f ~ -1-1 (m/cy). Furthermore as showed in this pa- 
per, the change of Sun's moment of inertia may also cause 
the significant contribution to orbital motion of planets. 
These facts may indicate the necessity of reconsideration 
of the astronomical system of units since AU is not only 
the conversion constant of two length units but also the 
value which characterizes GM of the Sun via the following 
relation, 

GMsun = fc'AUVd^ (24) 

where k = 0.01720209895 is the Gaussian gravitational 
constant and d means a day so that d = 86400 (s). 
Until now, several authors discussed the redefinition 
of the system of astronomical units e.g., Huang et al. 
(1995); Standish (2005); Klioner (2008); Capitaine and 
Guinot (2009). By using Eqs. (1) and (24), we obtain 
GMsun in SI units and so far GMsun is regarded con- 
ventionally as the "fixed value" in SI units. However the 
discussions in Noerdlinger (2008) and present paper imply 
that the effects of solar mass loss and the rotational an- 
gular momentum transfer of the Sun due to the radiative 
mass loss cannot be disregarded in terms of the growing 
observational accuracy and the definition of system of as- 
tronomical units. 

In present paper, we proposed the possible mechanism 
for explaining the secular increase of AU, nonetheless we 



need to verify the validity of our model by means of the 
some tidal dissipation models of the Sun. Moreover be- 
cause the existence of dAU / dt is not confirmed robustly in 
terms of the independent analysis of observation by other 
ephemerides groups, it is important not only to perform 
the theoretical researches but also to re-analyze the data 
and to obtain more accurate value of dA\J/dt adding new 
observations e.g.. Mars Reconnaissance Orbiter, Phoenix 
and forthcoming MESSENGER which is cruising to the 
Mercury. It also seems to be meaningful to use the obser- 
vations of outer planets as well, such as Cassini, Pioneer 
10/11, Voyager 1/2, and New Horizons for Pluto since it 
is more natural situation that the variation of moment of 
inertia of the Sun causes the orbital changes not only of 
inner planets but also of outer ones. Therefore in order to 
reveal the origin of secular increase of AU, it is essential 
to investigate these subjects in detail. 
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